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Abstract 

CN ■ We study a new class of kinetic equations describing nonequilibrium macroscopic 

^ . dynamics of soliton gases with elastic collisions. These equations represent nonlinear 

\^ ' integro-differential systems and have a novel structure, which we investigate by study- 

ing in detail the class of A^-component 'cold-gas' hydrodynamic reductions. We prove 
'""I . that these reductions represent integrable linearly degenerate hydrodynamic type sys- 

I tems for arbitrary N which is a strong indication to integrability of the full kinetic 

' equation. We derive explicit representations for the Riemann invariants and character- 

■ istic velocities of the hydrodynamic reductions in terms of the soliton gas component 

^ ! densities and construct a number of exact solutions having special properties (quasi- 

periodic, self-similar). Hydrodynamic symmetries are then derived and investigated. 



1 Introduction 



In the recent paper [38j V.E. Zakharov has put forward a programme for the construction 
of the theory of wave turbulence in integrable systems. Such an unconventional union of 
integrability and stochasticity has a clear physical motivation: nonlinear dispersive waves, 
while often being successfully modeled by integrable systems, could demonstrate very com- 
plex behaviour calling for a statistical description characteristic of the classical turbulence 
theories. 

One of the important problems arising in this connection is the description of "soliton 
gases" — random distributions of solitons which can be mathematically defined in terms 
of generalized reflectionless potentials with shift invariant probability measure on them (see 
e.g. [21] )• Due to isospectrality of the "primitive" microscopic evolution, the macroscopic 
dynamics of a homogeneous soliton gas is trivial (for strongly integrable systems, such as 
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the Korteweg - de Vries (KdV), nonlinear Schrodinger (NLS) or Kadomtsev-Petviashvili 
(KP-II) equations — see [3B]), namely, all the statistical characteristics can be specified 
arbitrarily at the initial moment and remain unchanged in time. However, if the soliton gas 
is spatially inhomogeneous, i.e. if the probability distribution function depends on the space 
coordinate, then nontrivial macroscopic dynamics occurs due to phase shifts of individual 
solitons in their collisions with each other. The approximate kinetic equation describing 
spatial evolution of the soliton distribution function in a rarefied gas of the KdV solitons, 
when these phase shifts can be taken into account explicitly, was derived by Zakharov back 
in 1971 [35]. 

Generalization of Zakharov's kinetic equation to the case of soliton gas of finite density has 
been made possible rather recently [5] and required consideration of the thermodynamic limit 
of the Whitham modulation equations associated with finite-gap potentials (note that the 
idea to use the Whitham equations for the modeling of a turbulent motion - a "deterministic 
analoque of turbulence" - was first proposed by P.D. Lax [25]). In the thermodynamic limit, 
the nonlinear interacting wave modes transform into randomly distributed localised states 
(solitons) and the modulation system assumes the form of a nonlinear kinetic equation. 
This new kinetic equation was extended, using physical reasoning, in [7] to other integrable 
systems with two-particle elastic interactions of solitons (i.e. when multi-particle effects are 
absent). 

The kinetic equation for solitons in general form represents a nonlinear integro-differential 

system 

ft + {sf), = 0, 



oo 
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Here /(?]) = f{ri,x,t) is the distribution function and s{ri) = s{ri,x,t) is the associated 
transport velocity. The (given) functions S{ri) and G{rj, n) do not depend on x and t. The 
function Girj^jj) is assumed to be symmetric, i.e. G{ri,fj,) = G{fi,r]). The choice 



S{v) = , G(r/,/x) = log 



7] + fi 



(2) 



corresponds to the KdV case, ipt — 6(p(Px + '^xxx = (see [5]). In the KdV context, f] > 
is a real-valued spectral parameter (to be precise, before the passage to the continuum limit 
one has Afc = —rjl, where Xk, k = 1, . . . ,N are the discrete eigenvalues of the Schrodinger 
operator), thus the function f{r],x,t) is the distribution function of solitons over spectrum 
so that K = f{j])drj = 0{1) is the spatial density of solitons. If k ^ 1, the first order 
approximation of ([T]), ([2]) yields Zakharov's kinetic equation for a dilute gas of KdV solitons 

The quantity S{ri) in ([1]) has a natural meaning of the velocity of an isolated (free) soliton 
with the spectral parameter r] and the function ^G{ri,fi) is the expression for a phase shift 
of this soliton occurring after its collision with another soliton having the spectral parameter 
fi < 7]. Then s{r],x,t) acquires the meaning of the self-consistently defined mean local 
velocity of solitons with the spectral parameter close to rj (see [7]). 

Theory of nonlocal kinetic equations of the form ([1]) is not developed yet. Possible 
approaches to their treatment were discussed in [2] in connection with special classes of exact 
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solutions for the Boltzmann kinetic equation for Maxwellian particles. The derivation of ([T]), 
([2]) as a certain (albeit singular) limit of the integrable KdV-Whitham system suggests that 
this new kinetic equation is also an integrable system, at least for special choices of functions 
G{ri,n). A natural question arising in this connection is: what is the exact meaning of 
integrability for the equations of the type ([1])? 

Integrability of kinetic equations has been the subject of intensive studies in recent 
decades. For instance, integrability of the collisionless Boltzmann equation (which is some- 
times called the Vlasov equation) can be defined in terms of two other closely connected 
(even equivalent in some sense) objects: the Benney hydrodynamic chain [3], [36], |I5] and 
the dispersionless limit of the Kadomtsev-Petviashvili equation ([221 [23]. It turns out that 
all these three different nonlinear partial differential equations possess the same infinite 
set of A^-component hydrodynamic reductions parametrized by N arbitrary functions of a 
single variable [El [IT] (we note that the solutions to these A^-component reductions are 
parametrized, in their turn, by another A^ arbitrary functions of a single variable). This 
property was used in Ferapontov & Khusnutdinova [TOl [II](see also [37], [12], [18], [27] ) 
when introducing the integrability criterion for a wide class of kinetic equations, correspond- 
ing hydrodynamic chains and 2+1 quasilinear equations. Moreover, it was proved in [29] that 
the existence of at least one A^-component hydrodynamic reduction written in the so-called 
symmetric form is sufficient for integrability in the sense of [TU]. Another possible approach 
to analyse an integrable kinetic equation is to use the fact that it possesses infinitely many 
particular solutions determined by the corresponding hydrodynamic reductions (see [27] for 
details). 

The distinctive feature of the kinetic equation ([T]) is its nonlocal structure, which repre- 
sents an obstacle to the direct application to it of the approaches developed in [29] and [27] . 
For instance, the possibility of an explicit construction of symmetric hydrodynamic reduc- 
tions (and even the existence of such reductions) for ([T|) are open questions at the moment. 
In this paper, we study a particular, yet probably the most important from the viewpoint of 
capturing the essential properties of the full equation, family of the 'cold-gas' A^-component 
hydrodynamic reduction of ([ID obtained via the delta-functional ansatz for the distribution 
function f(r], x, t) = J2iLi f\^^ i)^{v~Vi)- We prove that these reductions represent linearly 
degenerate semi-Hamiltonian (integrable) systems of hydrodynamic type (see [28] and [9]) 
and can be explicitly represented in the Riemann invariant form for arbitrary A^. The cor- 
responding characteristic velocities, conservation law densities and symmetries (commuting 
fiows) are described in terms of the so-called Stackel matrices determined by A^(A^ — 1) linear 
functions of a single variable. 

Integrability, for arbitrary A^, of the class of the hydrodynamic reductions studied in this 
paper is a serious argument in favour of integrability of the full nonlocal kinetic equation 
([I|), at least for certain choices of the functions S{X) and G{X,fi) in the integral closure 
equation. Of course, such an outcome does not look surprising for the particular choice 
([2D of S{X) and G{X,fi) corresponding to the thermodynamic limit of the integrable KdV- 
Whitham equations but our analysis suggests that the general integro-differential kinetic 
equation ([I|) is a representative of a whole new unexplored class of integrable equations with 
potentially important physical applications. 

The structure of the paper is as follows. In Section 2 we outline the derivation of the 
kinetic equations for the gas of the KdV solitons following the thermodynamic limit proce- 
dure of [5] and extending it to the entire KdV-Whitham hierarchy. We then introduce the 
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generalized equation ([T]), and in Section 3 consider its A^-component 'cold-gas' hydrodynamic 
reductions having the form of hydrodynamic conservation laws. We then formulate our main 
Theorem 3.1 stating that the hydrodynamic reductions under study are linearly degenerate 
and integrable (in Tsarev's generalised hodograph sense) hydrodynamic type systems for 
any N. Section 4 is devoted to the account of the main results of the theory of linearly 
degenerate hydrodynamic type systems. In Section 5 we prove the statement of the main 
Theorem 3.1 for the case = 3 by explicitly constructing the corresponding Stackel matrix 
and presenting expressions for the Riemann invariants and characteristic velocities in terms 
of the conserving component densities. We also construct two distinguished families of exact 
solutions (self-similar and quasi-periodic) to the 3-component reduction. In Section 6, the 
existence of the Riemann invariant parametrization of the cold-gas hydrodynamic reduction, 
via a single Stackel matrix, is proved for arbitrary N, which enables us to complete the proof 
of the main Theorem 3.1 for a general case. In Section 7, we derive explicit expressions for 
the Riemann invariants and characteristic velocities in terms of the component densities. 
And at last, in Section 8 we derive hydrodynamic symmetries (commuting flows) of the 
A^-component hydrodynamic reductions under study and then extract the family of linearly 
degenerate commuting flows. 

2 Kinetic equation for a soliton gas as the thermody- 
namic hmit of the Whitham modulation system 

We start with an outline of the derivation of the kinetic equation for the gas of the KdV 
solitons as the thermodynamic limit of the KdV- Whitham system following We then 
naturally extend this derivation to the whole Whitham-KdV hierarchy. 

Let us consider the Whitham modulation system associated with the A^-gap potentials 
UN{x,t) of the KdV equation. This is most conveniently represented as a single generation 
equation in the [T3] form: 

{dpN)t = {dqN)x , (3) 

where dpN and dqN are the quasimomentum and quasienergy differentials deflned on the 
two-sheeted hyperelliptic Riemann surface of genus A^ : 

2Af+l 

r : ;x2(A) = JJ (A - Xj) , A e C, A,- e R . (4) 

Ai < A2 < ■ ■ • < \2N < X2N+I 5 

with cuts along spectral bands [Ai,A2], • • • [A2j-i, A2j], [A27V+1, c>o]. We introduce the 
canonical system of cycles on F as follows (see Fig. 1): the a^-cycle surrounds the j-th cut 
clockwise on the upper sheet, and the f3j- cycle is canonically conjugated to a/s such that 
the closed contour f3j starts at A2j , goes to +00 on the upper sheet and returns to A2j on 
the lower sheet. 

The meromorphic differentials dpN and dq^ are uniquely deflned by their asymptotic 
behaviour near X = —00: 

- A > 1 : dpN^ - . t W2 ' ~ {-Xy/^dX (5) 
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Figure 1: The canonical system of cycles on the hyperelliptic Riemann surface of genus N . 



and the normalization 

/-^ 2N+1 

dqN = 0, i = l,...,N; CN = --^Xj. (6) 
3. i=i 

The integrals of dp^ and dq]^ over the a - cycles give the components of the wave number 
and the frequency vectors respectively 



rfpiv(A) = kj{\i, X2N+1) , j> dqN{X) = Wj(Ai, . . . , A27V+1) , j = 1,...,N . 

(7) 

Let Ai = —1, A2Af+i = 0. Following Venakides [31] we introduce a lattice of points 

1 ^ r^i > r/2 > • • • > ^iv ~ , (8) 

where 

-r7| = ^(A2,-i + A2,) (9) 



2 

are the centres of bands. 

We now assume that the spectral bands are distributed such that one can introduce two 
positive continuous functions on [0, 1]: 

1. The normalized density of bands if{ri): 

, , , number of lattice points in (77, ri + drj) 
if{ri)dri ^ — . 

That is, 



N{Vj - Vj+i) N J 



(10) 



2. The normalized logarithmic band width 7(77): 



liVj) = ^ogSj + O(^) , 6j = - \2j-1 ■ (11) 

The functions ipirf) and 7(77) asymptotically define the local structure of the Riemann surface 
F (jl]) for ^ 1. In other words, instead of 2N + 1 discrete parameters \j we have two 
continuous functions of rj on [0, 1] which do not depend on a;, t on the scale of the typical 
change of Aj's in ([3]), say Ax /S.t ^ I. 
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The existence of the continuous distributions ipirj) and 7(77) imphes the following band- 
gap scaling for N ^ 1: 

|band,| ~exp{-7(r/,)iV}, j = l,...,iV (12) 



The scaling (fT2|) has an important property: it preserves the finiteness of the integrated 
density of states as iV — )■ 00. The integrated density of states is defined in terms of the real 
part of the quasimomentum integral (see j2l]): 



A 

A/'jv(A) = ii?e j dpN{\') , A G [-1, 0] 



(13) 



Now, using ([7]) one can readily see that 

M 

27r 



1 

^N{\) = —y^k, if AG [A2M,A2M+l], M=l,...,iV, (14) 

/IT f * 

which is a particular (finite-gap) case of the general gap-labeling theorem for quasi-periodic 
potentials It is not difficult to show that the scaling (IT^ implies that kj ~ 1/A^ so the 
total density of states 

1 ^ 

i=i 

remains finite in the limit as — 00. For this reason we shall call the limit as — )■ 00, 
defined on the spectral scaling ( 1T2|) . the thermodynamic limit. 

We shall not be concerned here with the existence and the exact meaning of the thermo- 
dynamic limit for the finite-gap potentials UN{x,t) (which is a separate important problem) 
but shall rather directly consider this limit for the associated Whitham system ([3]). It is 
however, instructive to note that it follows from (1121) that in the thermodynamic limit the 
band/gap ratio vanishes for each oscillating mode (i.e. kj — Vj = 1,2. ...,A^), so the 
thermodynamic limit of the sequence of finite-gap potentials associated with the spectral 
scaling f lT2|) is essentially an infinite-soliton limit. It was proposed in |6] that this limiting 
potential should be described in terms of ergodic random processes and can be viewed as a 
homogeneous soliton gas (or homogeneous soliton turbulence - depending on which of the 
two "identities" of a soliton is emphasized: the particle or the wave one). Then it is natural 
to assume that the same thermodynamic limit for the associated Whitham system should 
describe macroscopic evolution of the spatially inhomogeneous soliton gas. Indeed, as we 
shall see, the thermodynamic limit of the Whitham equations turns out to be consistent (in 
the small-density limit) with the kinetic equation for solitons derived by Zakharov [52] using 
the inverse scattering problem formalism. 

We first note that A/at (A) defined by f|T3l) is a monotone increasing positive function so 
(iA/iv(A) is a measure supported on the spectrum of the finite-gap potential un{x) [21j. Next 
we introduce a 'temporal' analogue of the density of states (|T3|l by the formula 



A 

Vjv(A) = ^Re j dqN{\') , A G [-1, 0] . (16) 
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Then integration of the generating modulation equation ([3]) on the real axis of A from —1 to 
-if G [-1,0] yields 

dtdMN{-f) = d^dVN{-f) , 1]. (17) 

Thus the finite-gap Whitham-KdV system can be regarded as the system governing the 
evolution of the spectral measure. 

Now we consider the thermodynamic limits of dAfN and dVN which we denote as 

(iA/AT — 7- 7r/(?7)(i?7 , dVN ^ —'^f{v)^{v)df] as N ^ oo, (18) 

where the limit is taken on the thermodynamic spectral scaling f|T2l) . Since 'Kf{ri)dri is the 
limiting spectral measure, the function /(//) has the natural meaning of the distribution 
function of the solitons over the spectrum (the meaning of the function s{i]) will become 
clear soon). The functions /(?]) and s{rj) were shown in [6], [5] to be expressed in terms of 
the ratio a{ri) = (p{ri)/'j{ri) of the lattice distribution functions ( ITOi) . ( |TT|) by certain integral 
equations, which are then combined into a single equation directly connecting /(r^) and s(?7) 

1 



s{ri) = Arf + 

7] 



j log 


r7 + /i 







f{^x)[s{r^)-s{^l)]d^i. (i9) 



We stress that in the continuum (thermodynamic) limit given by equations ( JT8l) . (|T9l) the 
explicit dependence of the density of states on the spectral branch points \j disappears. The 
only 'reminder' of the hyperelliptic Riemann surface F (jlj) is the kernel In \ri + ^\/\ri — /i| 
which arises as the continuum limit of the off-diagonal elements of the period matrix B of 
the Riemann theta-function BAr(a;,t|B) defining, via the Its-Matveev formula, the finite-gap 
potential (see [M] and [S]). 

Thus, integral equation flTIJl) can be viewed as a local (in the x,t - plane) relationship 
between the functions firf) and sirf) characterizing the soliton gas. Let / ^ 1 be the 
characteristic length at which the change of functions /(?;), s{ri) is small (of an order l/l <^ 
1). Next, in the spirit of modulation theory (see [33], [13]) we assume that on a larger 
spatiotemporal scale. Ax ^ I, /S.t ^ /, we have f(rj) = f{ri,x,t), s{ri) = s{ri,x,t) and 
postulate, using (ITS]) , that 

dtdAfN irdtfir], X, t)dr] , dj^dVN -iTdr,f{r], x, t)s{r], x, t)dri . (20) 

Then modulation equation flTTl) assumes the form of a conservation equation for /, 

/* + (s/)x = 0, (21) 

which is clearly an expression of the isospectrality of the KdV evolution. Since p{x,t) = 
Jq fdr] is the density of solitons the quantity s(?7,x,t) can naturally be interpreted as the 
velocity of the soliton gas (or, more specifically, the velocity of a 'trial' soliton with the 
spectral parameter A = —r]'^ - see [20j). One can see from (fT9|) that this velocity is different 
from the velocity Af of the free soliton with the same spectral parameter. This difference is 
obviously due to the collisions of the 'trial' ?7-soliton with other '//' - solitons in the soliton 
gas. Indeed, for small densities p = J fdrj ^ 1 one can consider the second term in (1191) as 
a small correction to the free-soliton velocity and obtain to the first order in p 



s{r]) ~ Arf H — 

7] 



j In 


r] + p 




7]- p 



f{p)[Ar,^-Ap^]dp (22) 
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which is Zakharov's expression for the average velocity of a trial soliton in a rarefied soliton 
gas, obtained by taking into account the change in the soliton position due to phase shifts 
in its pairwise collisions with other solitons [35] . 

Equations f l2T|) and f lT9|) thus provide a self- consistent kinetic description of the KdV 
soliton gas of a finite density. 

The outlined procedure of the thermodynamic limit can be readily extended to the entire 
Whitham-KdV hierarchy, 

idpN)t^ = idqi;^h, neN, (23) 

where tt, is a number of the "higher" Whitham-KdV equation in the hierarchy (the original 
modulation equation (j3]) corresponding to the KdV equation itself has the number n = 1) 
and tn is the corresponding "higher" time, so that {dpN)t„tm = {d,pN)t^tn n ^ m. The 

meromorphic differential dqf^'' is uniquely defined by its asymptotic behaviour near A = — oo, 

~ {-Xy-^'^dX , (24) 



and the normalization 



jdq^^\X)=0, J = l,...,iV (25) 



analogous to 

Now, applying to equation (l23l) the outlined above procedure of the thermodynamic 
limiting transition we obtain the same transport equation (!2T!) for the distribution function 

/*„ + (s„/). = 0, (26) 
and the integral closure equation for s„ assumes the form 

1 



1 log 









f{fJ')[sn{v) - Sn{f^)]dn, (27) 



where C„ are certain constants. We won't need specific expressions for them here. Moreover, 
since the characteristic speeds of the commuting KdV-Whitham flows, and, therefore, the 
corresponding transport velocities s„ in the thermodynamic limit equation (127|) . are defined 
up to a constant factor, hereafter one can assume C„ to be arbitrary constants. 

We note that equation fl27p differs from (IT^ only in the first term corresponding to the 
free-soliton velocity. Also note that the 'phase-shift' logarithmic kernel in the integral equa- 
tion f l27|) is the same for all n, which is not surprising as the whole finite-gap Whitham-KdV 
hierarchy fl23|) is associated with the same Riemann surface, i.e. with the same period matrix 
B responsible for the form of the integral kernel in the limit. 

Now it is only natural to consider a generalization of the derived kinetic equations (l26l) . 
( 1271) by introducing in ( 127|) an arbitrary function S{ri) instead of the free-soliton velocity 
term and a symmetric function G{r), /i) instead of the 'phase-shift kernel' in the integral 
term. As a result, we arrive at the generalized kinetic equation ([1]), which will be our main 
concern in the remainder of the paper. 
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3 'Cold-gas' hydrodynamic reductions 



We introduce in ([T]) u{rj,x,t) = rif{ri,x,t), v{ri,x,t) = —s{ri,x,t) and consider an N- 
component 'cold-gas' ansatz 

N 

U 

i=l 

which reduces ([1]) to a system of hydrodynamic conservation laws, 



J2Ax,mV-Vi), (28) 



ui={u'v%, ^ = l,...,iV, (29) 
where the velocities and the 'densities' satisfy algebraic relations 

v' = ^^ + Yl - v'), = efc.. (30) 

Here ^ 

^i = S{r]i), eik = G{r]i,r]k), v' = v\x,t) = s{r]i, x,t) (31) 

ViVk 

Without loss of generality one may assume that en = for all i. 

In a two-component case, the above algebraic system fl30l) can be easily resolved for u^'^ 
in terms of w^'^: 

u' = u' = 32 

ei2 — ei2 f — 

Substituting (1321) into (1291) we arrive at the 

Lemma 3.1 (El & Kamchatnov 2005 [7J): Hydrodynamic type system ([13 ), (3^) for 
N = 2 reduces to a diagonal form in the field variables and v"^: 



.1 „,2„,1 „,2 „,1„,2 



(33) 



Remarkably, the hydrodynamic type system (l33l) is linearly degenerate because its charac- 
teristic velocities do not depend on the corresponding Riemann invariants. 

It is worth noting that system (l33l) is equivalent to the ID Born-Infeld equation (Born 
& Infeld 1934) arising in nonlinear electromagnetic field theory (see [33], |1]) 

(1 + vDVyy - 2(^xVyVxy + (1 - vl)Vxx = 0. 

As any two-component hydrodynamic type system, (133!) is integrable (linearizable) via the 
classical hodograph transform. However, for any > 3 integrability of the original system 
(1291) . (1301) is no longer obvious. As a matter of fact, most A^-component hydrodynamic type 
systems are not integrable for A^ > 3. Also, it is even not clear whether A^-component system 
(|29|) . fl30|) is linearly degenerate. It might seem that this system is simple enough for one 
to be able to verify these properties by a direct computation, using general definitions of 
linear degeneracy and integrability for hydrodynamic type systems [28], [3ll|32] (also see the 
next section). To our surprise, even the simplest non-trivial case A^ = 3 turned out to be 
complicated enough to require computer algebra to get the confirmation of our hypothesis. 

The identification of the system fl29l) . (!30|) for A^ = 3 as an integrable linearly degenerate 
hydrodynamic system can be considered as a strong indication that both properties (linear 
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degeneracy and integrability) could hold true for this system for arbitrary A^. Thus we 
formulate our main 

Theorem 3.1 N -component reductions [29\) . ^3D^ of the generalized kinetic equation (QP 
are linearly degenerate integrable hydrodynamic type systems for any N. 

To prove this theorem, we take advantage of the well-developed theory of integrable 
linearly degenerate hydrodynamic type systems [2B], [S]. 



4 Linearly degenerate integrable hydrodynamic type 
systems: account of properties 

In this Section we give a brief account of the general theory of semi-Hamiltonian linearly 
degenerate hydrodynamic type systems. 
A hydrodynamic type system 

Ui = v]{V)Ui, I = l,2,...,iV,j = l,2,...,iV (34) 



is called semi-Hamiltonian (see [3T1[32]) if it 

(i) has mutually distinct eigenvalues A*(U) 

det |A(5i -t;i(U)| = 0; (35) 

(ii) admits invertible point transformations U^{r), such that this hydrodynamic type 
system can be written in diagonal form 

rl = V\r)ri, t = l,...,N. (36) 

The variables r^{\J) are called Riemann invariants and V^(r) = A'^(U(r)) - characteristic 
velocities. Each Riemann invariant is determined up to an arbitrary function of a single 
variable -Rj(r*). 

(iii) satisfies the identity 

= '^^^^ (37) 

for each three distinct characteristic velocities {dk = d/dr^). 

A semi-Hamiltonian hydrodynamic type system possesses infinitely many conservation 
laws parameterised by A^ arbitrary functions of a single variable. Its general local solution 
for d^r^ ^ 0, i = 1, . . . , N is given by the generalised hodograph formula 



x + V\r)t = W\r), (3^ 
where functions W\r) are found from the linear system of PDEs: 



]yi _ yyj yi _ yj 



t,j = l,...,N, i^j. (39) 



Thus, the semi-Hamiltonian property (!37|) implies integrability of diagonal hydrodynamic 
type system in the above generalised hodograph sense. 
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It is not difficult to verify that commuting hydrodynamic flows of (136|) are specifled by 
the hydrodynamic type system 

ri = W^{v)ri, j = l,...,iV, (40) 

where r is a new time (group parameter) so that {rl.)t = {rDr imphes fl39|) . 

A sub-class of linearly degenerate hydrodynamic type systems is determined by the prop- 
erty 

diV' = (41) 

for each index i. It means that each characteristic velocity does not depend on the corre- 
sponding Riemann invariant r*. 

Theorem 4.1 (Pavlov 1987 [2H]): If semi-Hamiltonian hydrodynamic type system / f5^) 
possesses conservation laws ^2^) with = f/*(r) and f*(U(r)) = V^ir) then this system is 
linearly degenerate. These conservation laws are parameterized by N arbitrary functions 
of a single variable. 

Proof: The semi-Hamiltonian property (i.e. integrability) is given by the condition ([3] 
We introduce, following Tsarev [32], the so-called Lame coefficients by 



d,lnH, = ^^^, ^^k. (42) 

Suppose that some semi-Hamiltonian hydrodynamic type system fl5B]) can be written in the 
conservative form f l2^ with f*(U(r)) = y*(r). In such a case 

dkU'-r^ = dkiU'V')-rl 

Since r{x,t) is an arbitrary solution of (!36|) we obtain equations 

V'-dkU' = dkiU'V'). (43) 

If k ^ i, then 

i.e. each of the conservation law densities is determined up to an arbitrary function of a 
single variable Pi{r') (cf. (02]) and ( 14^]) ). 

U' = H,-P,{r'). (45) 

If k = i, then it follows from (H3l) that diV^ = i.e the system is linearly degenerate. The 
theorem is proved. 

Remark: A subset {u''} of the conservation law densities {U''} satisfying a given system 
of conservation laws (e.g. ( 12^ . (I5U]) ) is selected by flxing the functions Pk (e.g. Pkij'^) = 1 
— see fl9T]) in Section 7). 

It is worth noting that the statement converse to Theorem 4.1 is not valid in general. 
Indeed, let us consider the two-component system of conservation laws, 

= {UW{U\ U')U = {U\\U\ U')),. (46) 
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Suppose this hydro dynamic type system is hnearly degenerate, then it can be written in 
Riemann invariants r^{U^, U^), r'^{U^, f/^) as follows: 

rl = V\r\ryi, r\ = V\r\ryi 

where V^^'^(r) = f ^'^(U(r)). Let us introduce new conservation law densities = f/^ + t/^ 
and U"^ = — . Then the system of conservation laws f l46p assumes an equivalent form 

where the characteristic velocities 

[/V + t/V 2 f/V-f/V 

V = V = 

no longer coincide with V^^(r^,r^) and V^^(r^,r^). 

The full theory of linearly degenerate semi-Hamiltonian hydrodynamic type systems was 
constructed by Ferapontov in ^ using the Stackel matrices 

/ 0i(ri) ... <Pj,{r^) \ 



lN-2(^N 



r 



lN-I/N 



r 



\ 



(47) 



where are arbitrary functions (it is clear that without loss of generality one can put 

(/)^~^(z) = z and the number of arbitrary function reduces to A^(A^ — 2)). Then characteristic 
velocities of such linearly degenerate hydrodynamic type systems are given by the formula 

det AP^ 



det A 



(1) 



(48) 



where A-^^ is the matrix A without kth row and ith column. The family of the conservation 
law densities W corresponding to the semi-Hamiltonian system ( 136|) . ( HHl) is determined by 

(cf. m) 



det A 



(1) 



where Pi{r'^), i 



det A 

N are arbitrary functions. 



(49) 



Corollary 4.1 The system of conservation laws [29\) is a semi-Hamiltonian linearly de- 
generate hydrodynamic type system if and only if the densities and velocities w*(u) admit 
representations = ?7*(r) and w*(U(r)) = V^ir), specified by ( (7^ , ( f^^ , via N functions 



The functions r''{x,t) then satisfy diagonal system (E3j. 
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Proposition 4.1 (Ferapontov 1991 [9]): Semi-Hamiltonian linearly degenerate hydrody- 
namic type system ^3^, has N — 2 nontrivial linearly degenerate commuting flows 

r^. = V(^)(rH, j = l,...,N, k = 3A,...,N, (50) 

whose characteristic velocities are determined as (cf. fj48l) ) 

det A^''^ 

^(1) = ^^^- (51) 
detA« 

Any characteristic velocity vector W(r) = {W^{r),W^{r), . . . ,W^{r)) specifying linearly 
degenerate hydrodynamic flow ( l40i) commuting with (136|) . ( l48l) . can be represented as a 
linear combination of the "basis" characteristic velocity vectors ^{k) fl5T]) (including "trivial" 
ones V(2) = V (see ( HHj) ) and V(i) = 1) with some constant coefficients bk- Thus, for any 
component there exists a decomposition 

N 



W' = Y,hViky (52) 



k=l 



Theorem 4.2 (Ferapontov 1991 [9]): General solution r{x,t) of the semi-Hamiltonian 
linearly degenerate system is parameterized by N arbitrary functions of one variable 
fk{r^) and is given in an implicit form by the algebraic system 



X 



(53) 



yjm^ M = 34 N 



(note the change of sign for t compared to p] due to a slightly different representation 
of the diagonal system (l36ll in this paper). We also note that formulae (l53ll represent 
an equivalent of the symmetric generalised hodograph solution ( l38l) for semi-Hamiltonian 
linearly degenerate hydrodynamic type systems. 

It is instructive to introduce, following Darboux [8J, the so-called rotation coefficients 

Afc = ^, ^^k, (54) 



where the Lame coefficients are defined by (1421) . Then expression ( l37jl for the semi- 
Hamiltonian property assumes the form 

di(3jk = (3ji(3ik, i^ j ^ k. (55) 

Using linear system fl39l) can be related to another linear system 

diHk = /3ikHi, i ^ k, (56) 
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via the so-called Combescure transformation (see [8]) 

= §. (57) 

In other words, one can show (see |32]) that the ratio of any two solutions to f l56|) satisfies 
system fl39|) for the characteristic velocities of the commuting flows fHOj) . We note that the 
particular solution Hi of f l56p corresponding to the characteristic velocities Vi of the original 
system fl5Bl) is expressed in terms of the Lame coefficient Hi as 

Hi = V'Hi. (58) 

Of course, general solution Hi of system 0561) . as well as general solution of the 
generalized hodograph equations (1391) . is parameterized by arbitrary functions of single 
variable. 

Theorem 4.3 (Pavlov 1987 p8]): The class of the semi-Hamiltonian linearly degenerate 
systems of hydrodynamic type is selected, in addition to 153]) . by the set of restrictions on 
the rotation and Lame coefficients 

dilnHi,= d,\n/3ji (59) 

for any index j ^ i. 

Proof : Let us consider the Lame coefficients for the linearly degenerate systems. Using 
(|12D, dUD we have 

d,V' = d,\nH,-{V^ -V'), i^j, d,V' = 0. 
The compatibility condition dildjV^) = dj{diV^) implies that 

didj In Hi = dj In Hi ■ di In Hj, i^ j. (60) 
Now one can see that the l.h.s. of (160|) can be written in the form 

(H \ H H 

-^Pjij = PijPji + -^^iPji - -^l^jAHi. (61) 

On the other hand, the r.h.s. of (HI]) is nothing but the product /3ijf3ji. Now (]5^ immediately 
follows from f ]60]) and f ]6T]) . The Theorem is proved. 

Now, suppose that the rotation coefficients fj5^ for some linearly degenerate hydrody- 
namic type system are given. Then restrictions fj59]) determine not only the Lame coefficients 
(112]) but also all other solutions of ( 156]) associated, via ( 157]) . with the characteristic veloci- 
ties ( 1^ . ( 15T]) of the complete set of linearly degenerate commuting flows. Indeed, one can 
see that equations ( 15^ . ( 159]) actually represent N systems ordinary differential equations so 
that each system contains differentiation with respect to only one Riemann invariant. Thus, 
the general solution Hi of system (15^ . (159]) is parameterized by arbitrary constants (see 
Proposition 4.1). 

Let us introduce particular solutions H^'^^ of system (^5^, dSH]) such that (see (jl5]) . 

m) 

^ik) = ~, k = 2,3,...,N, 
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where Hi = i/f \ Hi = Hf> (see (^). As a matter of fact, V^^^ = V\ V'^-^^ = 1. Then ([59]) 
can be written in a shghtly more general form, 

di\YlPji = di\YlH\''\ 

- for any k and j 7^ i. 

Thus, the full class of linearly degenerate semi-Hamiltonian hydrodynamic type systems 
is determined by conditions (159|) . (IMI) and (155|) . We note that system (1591) . (15^ and (!55|l is 
an overdetermined system in involution. Its integration leads to the aforementioned set of 
particular solutions of ( |56ll that can be parameterized via a Stackel matrix (see (H71) . ( HHll . 
(119]) and dSI])). 



5 = 3: explicit formulae 

We now apply the construction outlined in the previous Section to the first nontrivial (from 
the viewpoint of integrability) case = 3 of the hydrodynamic reduction ( l29i) . ( l30l) . To 
prove our main Theorem 3.1 for = 3 we make use of Corollary 3.1. 

Let us suppose that hydrodynamic system of conservation laws (1291) . (1301) is linearly 
degenerate and can be written in a diagonal form (136]1 . i.e. we suppose that there exists an 
invertible change of variables r-'(u) , j = 1,2,3, such that system (129]) assumes a diagonal 
form 

rl = V^ir)ri, J = 1,2,3, (62) 

where V^-'(r) = f-^(u(r)). 

We introduce the Stackel matrix (147|) . which for A^ = 3 can be written in the form 

Ei(ri) 52(r2) 53(r3) 
A = I Ai(ri) A2{r^) As{r^) | , (63) 
1 1 1 

where Ak{z) and Bk{z) are arbitrary functions. 

Now, by Corollary 3.1, if system f l29]) . f pO]) is linearly degenerate and semi-Hamiltonian 
then its diagonal representation flB^ must have characteristic velocities in the form fH5]) . i.e. 
for A^ = 3 we have 

, _ B.jr^) - B.jr^) , _ B,{r^) - B,{r') , _ B,{r^) - B,{r') 

A2(r2) - A3(r3)' A3(r3) - Ai(ri)' Ai(ri) - A2(r2) ' ^ ^ 

Then, using f H^ the corresponding conservation law densities m'^ are found in terms of 
Riemann invariants as 

(65) 

where Pj{r^) are arbitrary functions and the determinant of the Stackel matrix is given by 

det A = Ai(ri)[i?2(r2) - B,ir^)] + A2ir^)[B,ir^) - B^{r')] + A,ir^)[B,{r') - B^ir')]. (66) 

Substitution of flM)) - (!66!) into (!30!) yields expressions for the functions Ak{z), Bk{z), 
Pk{z),k = 1,2,3. 
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Before we present these expressions, we note that it follows from (IMj) . ( 166|) that functions 
Bk{z) are determined up to a constant shift which is then translated into a certain shift for 
the functions Pk{z). It turns out that this shift can be removed by the simplest change of 
the Riemann invariants, {r^ + constant) )■ (although the relationships between the shift 
constants for Bk-, Pk and r'^ are rather cumbersome) so that we eventually obtain 

Ai{r') = r\ B,{r') = Cir\ z = 1,2,3, (67) 

where 

>. ^3^12 — ^2^13 > ^1^23 — ^3^12 > ^1^23 ~ ^2^13 /^oN 

Cl = , C2 = , C3 = , (68) 

ei2 — ei3 ^23 ~ ^12 ^23 ~ ^13 

Pi = r + €23 , P2 = r + ei3 , P3 = r + ei2 • (69) 

ei2 — ^13 ^12 — £23 ^13 ~ ^23 

Direct verification shows that the diagonal system (1621) . ( l64l) . (1671) . ( l68l) is indeed equivalent, 
via f l65|) . fl66|) . fl69|) . to the original set of conservation laws fl29|) . fl30|) . where f'^(u(r)) = 

Thus, system fl29|) . fl30|) is consistent with formulae f l64|) . fl65l) defined by Stackel matrix 
( l63l) . Therefore, by Corollary 3.1, the three-component hydrodynamic reduction ( 129|) . ( I30l) 
is a linearly degenerate semi-Hamiltonian (i.e. integrable) hydrodynamic type system. 

Remark: As we have seen, the outlined construction has some additional inherent "de- 
grees of freedom", namely, three arbitrary constants due to non-uniqueness of the Stackel 
matrix specifying a given linearly degenerate semi-Hamiltonian system. The full set of arbi- 
trary constants removable by an appropriate change of the Riemann invariants will appear 
later in Section 5 where we shall consider iV-component hydrodynamic reductions with ar- 
bitrary > 3. 

Using fl64|) -f l69|) we obtain explicit expressions for the characteristic velocities V'^ and 
densities u'^ in terms of Riemann invariants, 

W . C^ZLi!!!, ^■^=<4^. V''=<^!l^, (70) 

2 3 3 1 1 2 

detA ' ' detA ' detA ' ^ ' 

where 

det A = (Cl - C2)rV2 ^ _ ^^)^2^3 ^ _ ^^)^3^i _ (72) 

We note that, unlike in the case N = 2, algebraic system (15U|1 cannot be resolved for u'^ in 
terms of for any odd (cf. corresponding formulae in Section 2), because determinant 
of the matrix A of linear system f pOj) 

Au = b, 

where Ai^ = eik{v^ — f*) and hi = — ^i, equals zero due to its skewsymmetry. For instance, 
for A^ = 3, the consistency condition of this linear system (i.e. the condition that the rank 
of augmented matrix equals 2) is given by the relation 



£23 



{v' - v'){^, - v') + euiv' - v'){^, - v') + e,,{v' - t;^)(6 - v') = 0. (73) 
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Direct substitution of = V^{r) (170|) into (173|) shows that it satisfies identically. 

Using (ITTi) . ([72]), (1^ . ( 1^ one can express the Riemann invariants in terms of the 
densities u'' explicitly, 



^1 



(ei2 - ei3)(ei2ei3W^ + ei2e23W^ + ei3e23^^ + £23) 
[(^3 - ei)ei2 + (6 - 6)ei3]wi - (6 - e3)(ei2w2 + eia^s + 1^ 



2 _ (^23 - ei2)(ei2ei3^^ + £12623^^ + £13623^^ + £13) 

" [(6-6)£23 + (6-6)£i2K-(6-ei)(£i2ni + e23«3 + l)' ^ ^ 



^3 



(ei3 - £23)(£i2£i3^^ + ei2£23M^ + ei3£23^^^ + £12) 

[(6 - ^3)613 + - 6)e23]«3 - (^1 - 6)(ei3^^1 + 623^2 + 1^ 



Direct substitution shows that expressions (17i|) and ( !70!) are consistent with original algebraic 
system (I5UI) where f-' = V'' (r(u)). 

It is instructive to look at what happens to the diagonal system fl62p when the density 
of one of the components in conservation laws fl29|) . say u^, vanishes. One can see from ([71]) 
that if = (this corresponds to vanishing of P3 in fl7Tl) ) then the Riemann invariant 
becomes a constant, 

4i - 42 

so that the equation for satisfies identically and system ( B2]) reduces to its 2-component 
counterpart ( 15^ for 

m') = V\r^{u\ M^ 0)) , u^) = V^{r\u\ M^ 0)) , 

as one should expect. Similar reductions occur for = and = 0, which lead to 
= Tq = constant and = Tq = constant respectively. As a matter of fact, any function 
W{r^) is also a Riemann invariant so one can choose a new set of Riemann invariants say 
W = — Tq so that = when = 0, which could be useful for applications. 

Now we consider some special families of solutions to the linearly degenerate system ( l62l) . 
701). 



a) Similarity solutions 

One can see that, owing to homogeneity of the characteristic velocities (170]) as functions 
of Riemann invariants, system ( 162]) admits similarity solutions of the form 

' ^-rff), 2 = 1,2,3, (75) 



t° \t 

where a is an arbitrary real number and the functions /*(t), satisfy the system of ordinary 
differential equations 

(r(l) + r)— + «r = 0, 2 = 1,2,3. (76) 
ar 

Here the functions V^{\) are obtained from ( ffOl) by replacing r* with It is not difficult to 
see that, due to the structure of the characteristic velocities, the case a = implies only a 
constant solution /* = Iq, where IqjI"^, /q are arbitrary constants. If a 7^ 0, the general solution 
of (176]) can be found in an implicit form using the generalised hodograph formulae ( 153]) . where 
for iV = 3 we substitute, according to ([63]), ([67]), 0^(0 = = CfcC, 'PliO = MO = 
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To obtain similarity solutions ( !75|) one should use in (153|) = 1'^%-, where /3 = 2 + 1/a 
and Cj, 2 = 1,2,3, are arbitrary nonzero constants. Then the requirement that the functions 
/* must depend on r = x/t alone leads to the algebraic system 

r = ciCi(/'r+c2C2(/'r+c3C3(/'r, 
-1 = ci(/ir+c2(/^r+c3(/^r, (77) 



= ci(/ir-i+c2(/2r-^+c3(/^r-\ 

where 7 = — l/a and we have also replaced Cj/7 1— )■ Cj. Direct substitution shows that 
solution defined ( 1771) indeed satisfies system ( 1761) . 

6j Quasiperiodic solutions 

Another interesting type of solutions arises when one introduces in ( 153|) (for = 3) 

7 

m=l 

where Ei < E2 < ■ ■ ■ < E-j are real constants. Then, according to ( 1671) . solution ( l53l) assumes 
the form 

^1 ^2 ^3 



= i^£=/f-£=/f^^, (80) 



which resembles the celebrated system for the multi-gap (here - three-gap) solutions of the 
KdV equation. Unlike ( 1781) - ( IHOl) . however, the three-gap KdV solutions correspond to the 
Stackel matrix §^ with the rows A^iO = ^, ^fc(0 = , ^ = 1, 2, 3 [9j. 

Proposition 5.1. For any constants Ci 7^ C2 7^ Cs 7^ there exists at least one set 
{El, . . . , Eq} such that the solution r*(x, t), i = 1,2,3 described by (7^ - (EOj is quasi- 
periodic in X and possibly in t. 

We present here a sketch of the proof. Availability of the solution in the form (I78|) - 
(!80l) implies the existence of separate dynamics of r-'-s with respect to x and t. Indeed, 
differentiating ( 1781) - ( IHOl) with respect to x for fixed t one readily obtains 

^ = (r^-r'=)^^^S3, ^,j,fc = l,2,3, i^j^k, (81) 

where 

n(ri, r2, ra) = (Ci - C2)rV2 ^ _ ^^^^2^3 ^ _ ^^^^3^1 ^ ^ ^g2) 
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- see (172|). 

Analogously, differentiating (ITS]) - (18 Op with respect to t for fixed x one obtains 



^ = (Or^' - Ckr') ^ n ' ^,J,A; = 1,2,3, ^ ^ A: . (83) 

One can see that the flows ( ISTl) and ( l83l) are consistent with the spatio-temporal dynamics 
(1361) . (ITOl) . We also note that equations ( IHTl) . ( I83l) resemble Dubrovin's equations for the 
auxiliary spectrum dynamics in the KdV finite-gap integration problem (see, for instance, 
[26]). 

Let us now suppose that 

G [El, E2], G [E,, E4], G ^e], (84) 

so that all a/ Rjir"^) are real. The above condition ( 184|) means that the point p = (r^, r^, r^) G 
lies within the rectangular box Kijj^ G with the vertices at {Ei, Ej, Ej.), i,j,k = 

Now, for any set of the constants (i, (2, Cs there exists at least one box Kij^k = K* G M^, 
which is not intersected by the cone n(r^,r^,r^) = 0. That is, inside K* the denominator 
n(ri,r2,r3) in flHTl) never vanishes. 

Assume now that the 'initial' values of r^,r^,r'^ for some x = Xq belong to K* . Then 
it follows from f lHTj) that, under the x-flow (t = const), the point p remains inside K* and 
undergoes "elastic" reflections at the faces of K* as x varies (note that, since 7^ for 
j 7^ fc, the factor (r-' — r'^) in ( IHTl) never vanishes so the reflections occur only at the faces 
of K*). Therefore, the motion is quasi-periodic with respect to x as long as conditions 
(18^ are satisfied. Indeed, the system (IHTj) possesses two integrals ( !79|) and ( IHOj) outside the 
"resonant" points, where 11 = 0, so it specifies a quasi-periodic motion on a 3-torus provided 
conditions fl84|) are satisfied. Of course, if conditions (18^ are not satisfied at x = Xq the 
solutions r*(x) may blow up and not be quasi-periodic. 

The proof of quasi-periodicity of the t-flow is similar, however, there is an additional 
requirement that the factor {Cj^^ — Ck'^^) (182|) should not vanish for all r G K* which 
might impose additional restrictions on the choice of Ei (that is for some {Ej} the motion 
can be quasi-periodic in x but not in t). 

We note that the quasi-periodicity of the x- and t-fiows can be proved directly from the 
solution ( 1781) - ( IHOl) . however the outlined proof using the dynamical systems arguments is 
qualitatively more transparent and more readily yields the "resonant" restrictions for x- and 
t-fiows. 

6 Integrability of A^-component hydrodynamic reduc- 
tions 

Now we prove our main Theorem 3.1 stating that the A^-component 'cold-gas' hydrodynamic 
reduction ( l29l) . ( l30l) represents a semi-Hamiltonian (i.e. integrable) linearly degenerate hy- 
drodynamic type system. For that, according to Corollary 3.1, it is sufficient to show that 
the conservation law densities and the transport velocities admit parametric represen- 
tations (H9|) and (148|) . = f/*(r) and f*(U(r)) = V*(r), via A^ functions in terms of the 
Stackel matrix ( H7|) . 
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We suppose that hydrodynamic type system (1291) . (130!) can be rewritten in a diagonal 
form fl36l) . and, moreover, the characteristic velocities V^{r) coincide with the expressions 

Now, substitution of (HHj) . ( H9|l into f l30|l leads to the algebraic system 



N 



J2 ^^k{-l)'Pk det A^f ) = det Af ^ - 6 det t = 1, . . . , N , (85) 

fc=i 

for Pk{r^) and (l)\{r^). Here the matrix A^^^"* is the matrix A with first two rows and zth 
and /cth columns deleted. In the derivation of f lS^ we have used the determinant Sylvester 
identity (see, for instance, Gantmaher 1959) 

. (12) _ det Aj^^ det Af ) - det Af ^ det A^ 
'^''^^^'^ " det A 

Expanding the determinants, 

N N 



det Af ) = [{-ir'<Pl det A^f ] , det Af ^ = [(-1)'=+Vl det A 



(12) 
ik 



k=l k=l 

we rewrite equations f l85l) as nonlinear systems for cj)^ and Pk, where k,n = 1, . . . , N, 

N 

Y,{-1)\</^1 - ^^<Pl + t^kPk) det Ajf =0, ^ = 1, . . . , iV. (86) 

k=l 

We recall that 0^"^ = r^ 0f = 1. 

One can now introduce matrices 5i obtained from the matrix A by deleting the first 
two rows and the i-th. column, and adding the first row with the elements 0^ — + eikPk- 
Thus, each matrix 5i has dimension (A^ — 1) x {N — 1). Then the above set of equations (186|1 
can be rewritten as 

det(5i = 0, z = l,...,iV, (87) 
which implies that the rows of each of the matrices 5^ must be linearly dependent: 

N-2 

n=3 ^ ' 

k = l,...N, i = l,...,N -1, k^i, 

where Ci^k ai'e arbitrary constants. These conditions can be considered as N linear systems, 

for fixed k each. Since all these systems are consistent the functions 0^ and Pk can be found 

by solving system fl88l) . 

Constants Cj^i cannot be equal to zero since in that case, according to fHHj) . the velocities 
would become undetermined. Therefore, without loss of generality we can set Cj^i = 1 

and the number of free constants becomes A^(A^ — 2). Thus, the following Proposition is 

valid: 
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Proposition 6.1: General solution of system ^8M) is determined by solutions 



det Bir'' + det Bl 
det Bk 



Pk 



det Bl; 
det Bk 



(89) 



of N linear systems ^8^, where B^, B\ 



Bl and B\: 



are matrices with elements 



b] 




for / = 2 
for / = 3 if 
for / > 3 



ly^i + l 



(90) 



and h\ 



'ki — "fcl — "fcl — 



^mk , O/^- 



i+1 m 



i+1 m 
ki 




where Cm,i are arbitrary constants such that det Bk 7^ 0. 

Remark: The set of constants Ci^m for which det Bk = has Lebesque measure zero or 
requires a very special choice of the parameters rjk. The exceptional case is the following: 
the vectors ^, 1 and are linearly dependent which yields, according to the definition (|3T|) . 
a set of equations for the special values rjk- 

Thus, we have proved that all elements of the Stackel matrix (1471) depend linearly on 
Riemann invariants and these elements are determined from the algebraic system (!30|) up to 
N{N — 2) arbitrary constants removable by an appropriate change of the Riemann invariants 
(for instance, by a shift in the case N = 3). By Corollary 4.1, the existence of such a Stackel 
matrix automatically proves the semi-Hamiltonian and linearly-degenerate properties of the 
hydrodynamic reductions fl29l) . fl30l) . 

Now, our main Theorem 3.1 is proved. 

7 Riemann invariants and characteristic velocities: ex- 
plicit construction 

The construction described in Sections 3 and 6 provides a proof of the existence of Riemann 
invariants for system fl2I?]) . fl5U]) for arbitrary A^. The Riemann invariants are found to 
parameterize system fl2^ . fl5Ul) via the sole Stackel matrix, which, by Corollary 4.1, implies 
linear degeneracy and integrability of this system. Explicit representations for conservation 
law densities and transport velocities f * in terms of the Riemann invariants are given by 
Ferapontov [9] formulae (H9|) . ( HHj) where the entries 0^ of the Stackel matrix (147|) and the 
functions Pkij^) are defined by formulae ( 189|1 - (!90|) . Using the functions (pl one also obtains 
the generalised hodograph solutions 

The outlined procedure, while providing general theoretical framework for the study of the 
'cold-gas' reductions of the kinetic equation for a soliton gas, seems to be not very convenient 
from the viewpoint of practical calculations. It also involves N{N — 2) intermediate constants 
Ci^rn, which introduce an additional unnecessary complication. It is, thus, desirable to have 
more direct representations for the Riemann invariants and characteristic velocities, which 
will also be free from these intermediate arbitrary constants. 

We shall make use of the Theorem 3.1 and show that, once the linear degeneracy and 
integrability properties of system ( 129|1 . (130|1 are established, explicit relations between the 
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Riemann invariants r and the conserved densities u can be found by a straightforward 
calculation. The calculation will involve the properties of the Lame coefficients outlined in 
Section 4. 

First, without loss of generality we choose the following normalization (see ( H5|) ) 

= Hk, (91) 

where H^s are the Lame coefficients (l42l) . Now, using Theorem 3.1 we assume that hydro- 
dynamic type system ( |29l) . ( I30l) can be rewritten in a diagonal form ( l36l) . so that = f/*(r) 
and v*(U(r)) = V^{r). For convenience, in what follows we shall use small m's and v^s only, 
assuming that Uj = Uj{r) = U^{r), Vj = Vj{r) = Vj(r). 

To obtain explicit formulae for the Riemann invariants of the hydrodynamic reduction 
fl29|) . f l30|) we need first to prove its so-called "Egorov" property. 

Definition 7.1 (Pavlov & Tsarev 2003 [30j): Semi-Hamiltonian hydrodynamic type sys- 
tem ( l34l) is called the Egorov, if a sole pair of conservation laws 

dta{u) = d,b{u), dtb{u) = d,c{u) (92) 

exists. In this case (see ( f^^ and ^5^), 

dia = Hl dib = HiHi, diC = Hl (93) 

while corresponding rotation coefficients ( [5^ become symmetric, i.e. 

Lemma 7.1: Hydrodynamic reductions (!29|) . (!30l) are Egorov. 
Proof: We consider the sum of conservation laws (129|) . (l30l) 



(a + 5^efcmw'"(t^™-t^') 

.k=l \ mj^k 



(94) 



One can see that, since the matrix is symmetric, the last term in r.h.s. of flMjl vanishes. 
Thus, simplifies to the form 

However, the flux S^^m^ of conservation law (195!) is nothing but the density of another 
conservation law which can be obtained by the same summation but with the special weights 
6, i-e. 

Comparison with definition fl^?2]) implies that in our case 

which completes the proof. 

Now we formulate the following 



22 



Theorem 7.1: The Riemann invariants of N- component hydrodynamic reductions [2§\) . 
'3D\j can he found explicitly as 



t = l N. 



(97) 



Proof: 

For the sake of completeness of our construction we first show that the hnear degener- 
acy property f HTj) of system f l29|) . f l30|) readily follows from the (assumed) existence of the 
Riemann invariants r^. Indeed, differentiating fl5U]) with respect to the Riemann invariant 
r* and taking into account that (see f H2]) . f PT]) ) 



dAnu'' 



V — V 



we obtain the expression 



which reduces, on using (j44|) . to the form 



0. 



(9^ 



Equation f l98|) can only be satisfied if diV^ = for all i (otherwise the field variables in 
the algebraic system f l30|l would cease to be independent). Thus system ( 129|) . f pOj) is indeed 
linearly degenerate. 

Now, differentiation of algebraic system f l30p with respect to the Riemann invariant r'^ 
yields 



which reduces, with an account of f lMj) and the linear degeneracy property, to 



[V — V 



1 + ^ei^M™ j dk Inn* - ^eirndku" 



Since all characteristic velocities v'' are distinct, the expression in square brackets must 
vanish for any pair of indices i and k, i.e. we have 



dk In 



1 + E ^imU'- 



k^i. 



(99) 
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Integration of (p9|) yields 

^e,„M'" + i?,(rV = -l, (100) 

where Ri{r''), i = 1, . . . , N are arbitrary functions. 

We now differentiate f llOOp with respect to the Riemann invariants r* and r^, which gives, 
on using (!9T!) and (!54|) . 

J]e,^Ar« + i^UrO + i?.(r09, In i?, = (101) 

and 

J] + Ri{r')(3ki + eiudk In 5^ = (102) 

respectively. Substitution of (!96|) into (p3!) gives 

= J^Am + 9, In Hi, Hi = iiHi + - ii)Pim. (103) 

By expressing (9i In Hi from the above first equation, fllOll) and (11021) reduce to the form 
R^{r')H, = Ri{r')J2P^m - J^e.^Pim - R-{r'), 

(104) 

Substitution of the expressions Ri{r^)Hi and eimH^ into (llOOp yields to a set of constraints 
R'i{r^) = 1, i.e. Ri{r^) = r* + Oj, where are arbitrary constants. Since any function of the 
Riemann invariant is a Riemann invariant as well one can put without loss of generality that 
Ri{f^) = f''- Then (11001) reduces to fp7|) . The Theorem is proved. 

Taking into account Ri{r^) = r* and eliminating Hi from (I104p we arrive at the linear 
algebraic system 

{r'tkm - ^ik^im)Pim + {t't^ " eife)/3ife = Cifc, i 7^ k (105) 

m^i,k 

for the rotation coefficients while f ll04p reduces (cf. the first formula in f ll03p ) to 

*. = E(i-i?)/'--7' w 

Now, using f l58|) . fl9T]) . f ll03p and (I106p we formulate the main result of this Section: 

Algebraic relations f pOj) can &e resolved in a parametric form: 



u 



E (i - 1?) ft". - "- = 6+Ek.»-«')^. (107) 
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where the symmetric rotation coefficients are found from the linear system (see (11051) ). 



XI ('^'^krn - eikeim)Pim + ij'r'' - e'^^^) f^ik = Uk, i ^ k. (108) 

m^i,k 

Thus, we have obtained unambiguous expressions for the Riemann invariants and character- 
istic velocities of the A^-component hydrodynamic reductions (|5U]) . 

In particular, for = 3 we have from (I107p . (llOSp explicit expressions for the rotation 
coefficients and characteristic velocities 

3 

p _ r ei2 — ei3C23 

^1^2^?, _ 7-1^2^ — ^"^^13 — ''"^£12 + 2612613623' 
2 

r 613 — ^12^23 

''^13 ~i 9 ^ 1 2 9 2 q 2 i i-i ' (109) 

^ij.z^a _ T-J^e^g — r"^6i3 ~ ''"^12 + ^612613623 

P23 '- 



j.1^2^3 - rl6^3 - r26f3 - T^ef^ + 2612613623' 



6 (e23 - r^r^) + 6 (^^^£12 - £13623) + "^3 (r^6i3 - 612623) 



6 


623 — r'^r^ + 613 (r^ 
(6^3 - rV^) + ^3 (rl623 - 


- £23) + ei2 {r^ - £23) 

- ei2ei3) + 6 ('"^eia - 


£13^23) 


^3 


6^3 — r^r^ + 613 (r'^ 
(6^2 - rV2) + ^2 (ri623 - 


- 613) + e23(r^ - 613) 

- eueis) + 6 (^^^£13 - 


ei2e23) 



(110) 



6^2 - + 6i3(r2 - 612) + 623(r^ - 612) 

One can observe that formulae ( II 101) do not coincide with representation (170|) . (!68|) obtained 
earlier for the same family of the characteristic velocities. The reason is that the two repre- 
sentations correspond to different choices of the Riemann invariants (we recall one more time 
that any function of a Riemann invariant is a Riemann invariant as well). The relationship 
between these two equivalent sets of the Riemann invariants is obtained by equating the 
characteristic velocities (I70|) and (IllOp . As a result we get 





(ei3 


- ei2){e23r^ - 612613) 




(6 




+ (6 - 6)ei2 + (^1 - 


6)ei3 




(ei2 


- £23) (£13?"^ - ei2e23) 




(^3 




+ (6 - 6) £23 + (6 - 


6)ei2 




(e23 


- ei3)(6i2r^ - ei3e23) 






- ^2^ 


+ (6 - ^3)613 + (^3 - 


6) £23 



(111) 



Here by r^, r^, we denote the 'old' Riemann invariants as in ( [70|) . 



8 Commuting hydrodynamic flows 
8.1 General explicit representation 

Commuting flows to semi-Hamiltonian linearly degenerate system (!29|) . (!30|) are deflned 
in terms of the Riemann invariants by equations (HOjl . (139|1 . We recall that, according to 
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Proposition 4.1, only — 2 of the commuting flows are linearly degenerate (excluding the 
'trivial' flows specifled by linear combinations of the constant characteristic velocity 1 and the 
characteristic velocity v of the original flow 0361) ). The general solution of the generalised 
hodograph equations fl39l) specifying commuting flows was obtained by Ferapontov |9] in 
terms of the Stackel matrix entries (see Theorem 4.2). Here we are interested in a more 
explicit representation of the commuting flows for the speciflc system fl29|) . fl30|) . For that, 
instead of integrating system ( l39l) . we take advantage of the fact that our linearly degenerate 
system ( l29l) . ( l30i) is Egorov. In that case, the commuting flows can be found explicitly. 

We flrst observe that any conservation law density h for linearly degenerate hydrodynamic 
type system (129!) can be represented in the form (see (145|) or (149|) ) 



N 

h = Y,n''Pk{r^), (112) 

fc=i 

with arbitrary functions Pki^^) of a single variable. Then we make use of 

Lemma 8.1 (Pavlov & Tsarev 2003 [30]): All commuting flows (2^? ^he Egorov 

case are specified by the expression (see ( [57| j, ( fPgj) ) 

Hi diG 

Substituting (1112^ . ( l96l) into ( 1113^ and using the flrst formula from (1103^ we obtain an 
explicit representation for the characteristic velocities of the commuting flows (HOl) . 



W' = P,{r') + y- Pl{r') + YiP^ir^) - P,ir'))/3,^ . (114) 




We recall that Pfc(r'=), k = 1,...,N are arbitrary functions and the dependence of the 
rotation coefficients (3im on the Riemann invariants is found from system (11081) . 

If Pkij^) = 1, (I114p reduces to = 1; if Pkij^) = ^k, it reduces to the second formula 
in (11071) . i.e. to hydrodynamic reduction (!29|) . (!30l) itself. 



8.2 Linearly degenerate commuting flows 

To extract the family of linearly degenerate commuting flows from general representation 
(11141) we formulate the following 

Lemma 8.2: For the linearly degenerate commuting flows each function Pi{r^) in (TJ^ 
is linear with respect to the corresponding Riemann invariant r\ 

Proof: The condition diW^ = of linear degeneracy of the commuting flow implies, on 
using ([SSD and ([IDS]), that Pl'{r') = 0. 

We now consider the representation for the family of linearly degenerate commuting 
flows suggested by the form of the kinetic equations ([26]) . (I27j) for the KdV hierarchy. Im- 
portantly, the whole KdV kinetic hierarchy (l26l) . (1271) is characterised by a single integral 
kernel, G{ri, fi) = \ia\{ri — fi) / {rj + (which is consistent with the fact that all equations of 
the original flnite-gap Whitham hierarchy are associated with the same Riemann surface). 
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This suggests that there could exist a family of commuting flows to general nonlocal kinetic 
equation ([T]) having the form 



oo 



~s{v) = S{v) + - [ G{r^,f^)f{f,Mf^) - 5(r/)]ci/z, ^^^^^ 



V 



where S{r]) is an arbitrary function. Although verification of commutativity of the kinetic 
equations ([1]) and ( IllSp is beyond the scope of the present paper, it is clear that, if these 
equation do commute, this must be manifested on the level of hydrodynamic reductions as 
well. Having this in mind, we consider the iV-component hydrodynamic reductions to flllSp 
obtained by the familiar delta-functional ansatz (1281) and try to see if they commute with 
the original reductions f l29|) -f l3T|) . 

First we notice that equation flllSp is, essentially, the same kinetic equation ([1]) but with 
a different time variable and different "free soliton speed" function S{ri). Now, since we 
have proved integrability of the linearly degenerate hydrodynamic reductions fl2I?l) - flSTl) in a 
general form, we automatically have that analogous A^-component hydrodynamic reductions 
of flllSp must also be integrable linearly degenerate systems. It should be noted that, since 
the function S{ri) is arbitrary, the set {^i, . . . ,^Ar} of its values = S{rij) can be viewed 
as a set of N arbitrary numbers, and the corresponding 'cold-gas' hydrodynamic reduction 
becomes (cf. ([29]) -(EI])) 

< = (nV),, 2 = 1,..., AT, (116) 
where the velocities = P and the conservation law densties satisfy algebraic relations 

N 

v" = ^i + ^eiku''{v'' -v'), eik = eki, en = 0, (117) 

k=l 

and ejfc are the same as in flHTl) . 

According to Theorem 3.1, equations f lll6p . flll7p can be represented in the Riemann 
form 

rl = ij\r)rl, ^ = 1, 2, AT; /c = l,2,... (118) 

where the dependence {;*(r) of the characteristic velocities on the Riemann invariants is 
determined by the same formulae (I107P with the only difference that, one now replaces C,j 
with ^j, i.e. we have 

v' = i^ + J2(^m-^.)^. (119) 

Indeed, representation flll9p is a straightforward consequence of fll07p since the rotation 
coefficients /3im and Lame coefficients Hi do not depend on the parameters (see f llOSp . 
first formula in (11070 . and normalisation fl9T]) ). It not difficult to see that commutativity 
relationships (see fl39]) ) 

«,J = 1,2,3, t^j, (120) 



yi — — 

are identically satisfied. Thus, we have proved the following 
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Lemma 8.3: Linearly degenerate semi-Hamiltonian flows U16\) . (Ill) and [29^] . [W^ 
commute for any N . 



In conclusion we note that, although we have proved integrability of the 'cold-gas' hy- 
drodynamic reductions ( l29l) . ( l30l) for an arbitrary choice of the functions 5'(?7) and G(?7,/i) 
in the original kinetic equation ([1]), one can expect that integrability of the full equation 
([1]) would require some additional restrictions on the integral kernel Gij], fj,) (other than just 
symmetry) . 
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